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We report a recent lattice-QCD calculation of nucleon and pion electromagnetic form factors and 
'~J^ nucleon axial form factors, with special emphasis on large Q^. Conventional lattice form-factor 

calculations can only reach about 2.5 GeV^, but in this work the transfer momentum is pushed as 
large as 6 GeV^. Here, we demonstrate the results on 2+1 -flavor anisotropic clover lattices for the 
nucleon and pion, comparing with low-g' quantities, such as Dirac and Pauli radii, anomalous 
^""^ magnetic moments, gA and Ma- Our approach can be applied to isotropic lattices and lattices with 

smaller lattice spacing to achieve even larger-Q form factors. The form factors are processed 
to obtain transverse charge and magnetization densities across 2-dimensional impact-parameter 
space. These measurements could give important theoretical input to experiments, such as those 
of JLab's 12-GeV program, and provide insight into hadronic structure. 
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1. Introduction 

Nucleons are the building blocks of ordinary matter and pions are the lightest hadrons existing 
in our world. In the Standard Model, they are composed from fundamental particles, quarks and 
gluons; non-trivial interactions of quarks and gluons within the confined space of hadrons deter- 
mine the properties of these hadrons. Due to confinement, we cannot directly observe how the 
quarks and gluons compose the pion and nucleon; rather, experiments smash hadrons into differ- 
ent final states (as in deep inelastic scattering) to provide us with understanding through observ- 
ables, such as the parton distribution functions (PDFs). Processes such as electron elastic scattering 
with hadrons via photons give us electromagnetic form factors, which contain information in the 
momentum-space impact plane. Adding up information from various experiments, we can recon- 
struct these processes to gain a better understanding of how the quarks and gluons shape the nucleon 
and pion. 

Understanding form factors at high transfer momentum (Q^) is particularly interesting. It not 
only helps to reveal the nature of hadrons at short distances but also challenges models whose pa- 
rameters are calibrated in the \ow-Q^ regime. Experiments such as those in the 12-GeV upgrade at 
Jefferson Lab will soon provide higher precision and higher momentum-transfer form-factor data; 
progress will be made not only in the nucleon sector but also for pions. Theoretically, perturbative 
QCD (pQCD) should perform better at larger Q^, and we would like to get an understanding of 
when pQCD starts to become reliable. However, pQCD seems to fail to describe recent BaBar 
results for 7*7— >• Tl" over a wide range of transfer momenta 4 GeV^ < <40 GeV^ [1, 2], sug- 
gesting that there may be a large gap before pQCD can describe the experimental data correctly. 
A nonperturbative QCD approach would be able to provide the QCD theory results to connect the 
missing dots if one could use such an approach for sufficiently large Q^. 

The techniques of lattice QCD have been applied to such varied phenomena as the spec- 
troscopy of heavy-quark hadrons, the tower of excited baryon states, flavor physics involving the 
CKM matrix, hadron decay constants and baryon axial couplings. In many cases, it can provide 
higher-precision data from the Standard Model than what can be measured experimentally. By 
discretizing spacetime into a four-dimensional grid with a fixed lattice spacing and finite volume, 
we ai^e able to compute the path integral (in terms of discretized versions of the QCD Lagrangian 
and operators) directly via numerical integration, providing first-principles calculations of the con- 
sequences of QCD. The four-dimensional lattice size, lattice spacing, the quark masses in the sea 
sector and precision of the integration have been greatly improved due to the increasing compu- 
tational resources and more efficient algorithms over the past decade. Preliminary calculations 
at physical pion masses are being explored. Lattice calculations with sub-percent precision have 
been achieved for many flavor-related quantities and used to guide tests of Standard Model and 
to probe possible new physics. Calculations of hadron structures and interactions require more 
computational resources but should soon enter an era of precision. 

Hadron form factors have been calculated on the lattice by many groups, and calculations are 
still ongoing; see Refs. [3, 4, 5, 6] and references within for a few examples of recent review arti- 
cles. Recently, lattice calculations have also been used to calculate nucleon transition form factors 
involving excited nucleons [7] and experimentally challenging quantities, such as the hyperon ax- 
ial couplings [8], for the first time. However, the typical range in lattice calculations of hadron 
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form factors is less than 2.5 GeV^ on a typical lattice spacing a f« 0.12 fm. When one attempts 
higher-2^ calculations, they suffer from poor signal-to-noise ratio; for example, see a case study 
for the pion in Ref. [9]. 

There are two main obstacles in lattice QCD to achieving higher-2^ form factors. One limi- 
tation is the lattice spacing (a) available to be used in the calculation, since the momenta allowed 
on the lattice (assuming periodic spatial boundary conditions) are where L is the length of 

the spatial volume while n ranges 0, 1, V2, \/3, etc. Thus, the finer the lattice spacing, the higher 
momentum transfer one can achieve. However, very fine lattices may not be available in the next 
couple years and the spatial volume needs to be sufficiently large to avoid finite-volume effects, so 
we must instead increase the momentum in lattice units. In many existing calculations, the signal- 
to-noise ratio in form factors becomes quite poor before reaching high n at a given lattice spacing. 
We have examined the typical steps in a form- factor calculation to find the main causes of this noise 
and to develop solutions. The traditional method simplifies the complicated form-factor analysis 
by working in a regime where only the ground-state signal remains. This greatly simplifies the 
analysis, since taking ratios of the three-point correlator with current-insertion and the hadron two- 
point correlators will yield a constant in Euclidean time, which can be fit to a simple plateau. To 
suppress excited states in the limited time extent of the lattice box, we usually use interpolating op- 
erators that have broad extent in space, called "smeared"; this gives the same ground-state signal at 
eaiiier Euclidean time. The parameters of the smearing are usually tuned by looking at the hadron 
masses at rest or checking the form factors at zero or very low momentum transfer. However, as we 
increase the magnitude of the momentum transfer carried by the currents, the momentum earned 
by the initial or final hadron state must also increase (because of momentum conservation). The 
smearing parameters used to filter out the excited states will also filter out the boosted hadron sig- 
nal. Thus, in a system cai^efuUy tuned to include only the ground state, the signal will die out after 
a few momentum points when using fixed smearing parameters throughout the calculation. The 
solution to this problem is simple: we can vary the smearing parameters as the current momentum 
changes but keep a close eye on the contamination from excited states for each momentum; or we 
can take not only the ground state but also the excited state into the analysis, extracting the ground 
state using multiple smearing parameters or operators that has good overlap with various different 
states. 

We first investigated this idea using a quenched lattice with pion mass 720 MeV, calculating 
the J-quark connected contribution to nucleon form factors. For this demonstration, we fixed the 
final-state nucleon at rest so that the initial nucleon would carry the momentum transfer out at the 
inserted current. Ref. [10] shows the results from both methods: In the low-momentum region, we 
find the two approaches show consistent results with compatible statistical error bars. However, as 
soon as the momentum increases above the rest excited-state energy set by the smearing parameter, 
the initial nucleon operator contributes no further signal and presents a great deal of noise. At that 
point, some of our parameters that couple well to the rest excited state become dominant for the 
boosted nucleon ground state. These contribute to the form factors in the new method, and the 
signal-to-noise ratio remains relatively stable. Therefore, we can reach higher momentum transfer 
than the traditional fixed-operator, ground-state-only approach. 

In this proceeding, we will demonstrate the results of this improved method, focusing on re- 
sults from dynamical lattices. We will not only cover nucleon electromagnetic form factors, which 
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Figure 1: Isovector Dirac (left) and Pauli (right) form-factor extrapolation using Nf ^2+1 calculations in 
this work. The shaded bands show the simultaneous fit to the dependence of each pion mass points (875, 
580 and 450 MeV from top to bottom), where the color of the band matches the color of the points. The 
lowest gray band is the extrapolation to the physical pion mass. 

are more commonly calculated and better covered by experimental data over Q^, but also the less 
known large-Q^ region for nucleon axial and pion form factors. At the end of this proceeding, we 
will propose a step-scaling method for a future calculation of form factors at even higher momen- 
tum transfer, free of the systematics introduced in a calculation at fixed lattice spacing. 

2. Numerical Results 

The main body of our work uses the 2-1-1 -flavor anisotropic lattices generated by the Hadron 
Spectrum Collaboration (HSC) [11, 12]. (Note that the same techniques will work equally well 
for isotropic lattices.) These lattices use Symanzik-improved gauge action with tree-level tadpole- 
improved coefficients and anisotropic clover action [13]. The gauge links in the fermion action are 
3-dimensionally stout-link smeared with a small smearing weight p =0.14 and Up =2 iterations. 
The renormalized gauge and fermion anisotropics are around = 3.5 (that is, = 3.5at), and the 
inverse of the spatial lattice spacing is about 1.6 GeV. The quark propagators are calculated under 
antiperiodic boundary conditions in the time direction, while the spatial ones remain periodic. 
Quark propagators on the lattices are evaluated for source and sink operators with five Gaussian 
smearing parameters: a G {0.5, 1.5,2.5,3.5,4.5}. We use four final momenta (^/ = ^{0,0,0}a~^ 
^{-1,0,0}^"', ^{-l,-l,0}a"\ ^{-2,0,0}a"^) and vary the initial momentum over all pi = 
^{nx,ny,ny}a^^ with integer n^^,; nl + riy+n^ < 10 for current matrix elements. For more 
details on the parameters used in our work, please see the Ref. [10]. 

2.1 Nucleon Dirac and Pauli Form Factors 

From the HSC ensembles [11, 12], we use the 16^ x 128 lattices with pion masses of 875, 580 
and 450 MeV to calculate nucleon electromagnetic form factors. We extrapolate through a 
simple parametrization as used in interpolating experimental data: 
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where 1 = 0^/ (4M^). We simultaneously fit the m^j and (or t) dependence of the lattice data, 
expanding each 2^-fit parameter in terms of the pion mass: a,- = af'^ +a'pm\. The right-hand 
side of Fig. 1 shows an example of the extrapolations on the isovector Dirac form factor in the 
dynamical ensembles. The /Aof for F^j are 1.4 and 2.0, respectively. The lowest line/band 
represents the extrapolated form factors at the physical pion mass. We find good agreement between 
the Dirac form factor extrapolation at the physical pion mass and the experimental values; that is, 
the curvature (which is proportional to the product of squared nucleon mass and the mean-square 
charge radius) as a function of fits nicely with interpolating forms derived from experimental 
data. However, the Pauli form factor at physical pion mass suffers slightly larger ;^^/dof, mostly 
due to the failure to reproduce the (anomalous) magnetic moment, which would require smaller- 
momentum data ( < 500 MeV) that are not accessible with the current boundary conditions. At 
large momentum transfer, it becomes more consistent with the experimental form. 

To compare with other existing form-factor calculations, we can examine quantities in the 
small-2^ region, such as the isovector Dirac and Pauli mean-squared radii, from the isovector 
electric form factors FI2 via 



Most groups have studied radii with the dependence over ranges 0.5-2.0 GeV^ and found the 
extracted radii to be independent (within the statistical error bars) of choice [14, 15, 16]. The 
Dirac and Pauli mean-squared radii from the dynamical ensembles are summarized in the upper 
row of Fig. 2 along with other A'y = 2 + 1 lattice calculations and lowest-order heavy -baryon chiral 
perturbation theory (HBXPT) using experimental inputs [17]. A dimensionless quantities, Am^{r\) 
versus mj^L is also shown in Fig. 2, in case the discrepancy is caused by setting the lattice spac- 
ing. Our results are nicely in agreement with isotropic Nf = 2 + \ calculations having various 
sea and valence fermion actions; this demonstrates the universahty of the lattice-QCD calcula- 
tions. However, both mean-squared radii are a factor of 2 smaller than the experimental values, 
and HBXPT suggests a rapid increase of the radii for pion masses between 250 and 140 MeV; 
near-future calculations will be able to check this out. The (anomalous) magnetic moment, defined 
as Ktv = F2{Q^ = 0), is also shown in the right-most figure in Fig. 2. 

Since the elastic form factors contain information about the spatial structure of the nucleon, 
we can convert our data into a description of its densities in the impact plane. Having form factors 
at larger momentum transfer can greatly influence our understanding of the short-distance hadron 
structure. Due to the relativistic effects of the transferred momentum on the wavefunction of the 
nucleon, we cannot use a simple three-dimensional Fourier transformation without some recourse 
to models. Instead, we use the model-independent formulation of Ref. [19, 20] in terms of densi- 
ties in a two-dimensional plane transverse to an infinite-momentum boost. The transverse charge 
density (p) is defined as the Fourier transform of the form factor in such a plane: 




(2.2) 




(2.3) 
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Figure 2: Summary of the isovector Dirac (upper left) and PauU (upper right) mean-squared radii from all 
currently existing Nf = 2+1 nucleon electromagnetic form-factor calculations [14, 18, 15, 16, 10]. The 
dashed line indicates the leading-order HBXPT prediction, (lower left) Summary of the dimensionless 
products of Dirac radii and nucleon mass from each Nf = 2 + 1 lattice calculation, as a function of m^jL. 
The horizontal Une indicates the experimental values, (lower right) Summary of the normalized isovector 
anomalous magnetic moments from all currently existing A^^ = 2 + 1 nucleon electromagnetic form-factor 
calculations [14, 18, 15, 16, 10], including this work. 



where bold vectors b and q lie in the transverse plane. Equivalently, 

p{b) = j^^MbQ)F,{Q\ (2.4) 

for scalar b, where Jq is a Bessel function. We can perform this integral numerically, using the 
Pi{Q^) obtained by extrapolating our fit form to the physical pion mass. (Similarly, one can ex- 
amine the magnetic density (Pm) by integrating over J\ Bessel function and F2 form factor.) In 
Ref. [10], we demonstrate the importance of the high- 2^ form factors by restricting the data set 
to the region <2 GeV^ (which is the upper limit of transfer momentum for many lattice-QCD 
calculations) and comparing the result to the density obtained from using all the available Q^. The 
impact is significant in the central core, as we anticipated; omitting information about large transfer 
momenta results in a deviation in the density around 25%. (See Fig. 13 in Ref. [10].) Transverse 
charge density and magnetic density in infinite-momentum frame are shown in Figs. 14-16 in 
Ref. [10]. 
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Figure 3: The transverse charge densities in a polarized nucleon. The left- (right-) column figures present the 
results of proton (neutron) with b (lower) and b^^y (upper) in fm. In the upper-row figures, black indicates 
near-zero values, and purple, orange and white are increasingly positive. For the right figure, the same 
color-scale applies, slightly negative values are gold and the most negative are green. In the lower-row 
figures, the blue band indicates the densities from our lattice 2h-1 -flavor calculation, while the dashed line is 
a parametrization interpolating and extrapolating from available experimental data [21, 22]. (Note that the 
neutron experimental data is extrapolated to larger region). 

We can further look at the transverse charge densities in a polarized nucleon [23] via 

PT{b)=p+sin{^) r$^Ji{bQ)F2{Q^), (2.5) 
Jo IkMn 

where p{b) is obtained from Eq. 2.6. Results for the deuteron [24] and nucIeon-Roper transi- 
tion [25] have also been studied with experimental data. Figure 3 shows the results for the proton 
(left column) and neutron (right column) in one-dimension (lower row) and the two-dimensional 
impact plane (upper row) using our lattice inputs. Our lattice proton density is similar to the ex- 
perimental one, while the neutron has more significant deviations. This is likely due to the lack of 
disconnected diagrams in the form factors, which are closer to the magnitude of a typical neutron 
form factor but smaller than the proton one. Also experimentally, the neutron electromagnetic form 
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Figure 4: (upper left) Nucleon isovector axial form factors using 3 pion masses at 875, 580 and 450 MeV. 
The lowest gray band is the extrapolation to the physical pion mass. The dashed line is a dipole form using 
Ma = 1.03(2) GeV extending beyond the region of the available experimental data. gA (upper right) and 

^N\Jr\ (lower left) obtained from the same ensembles, and comparisons with previous Nf — 2+\ results, 
(lower right) Polarized distribution in a longitudinally polarized neutron. 



factors only known to 1.5 GeV^; thus a majority of the form-factor inputs are based on extrap- 
olation to larger region. Note that the asymmetry in the distribution for a polarized nucleon 
is due to the relativistic effect of boosting the magnetic moment of the baryon. This induces an 
electric dipole moment that shifts the charge distribution. 

2.2 Nucleon Axial Form Factors 

A similar approach can be applied to nucleon axial form factors, where experiments mostly de- 
rive results from neutron beta decay or pion form factors where various theoretical models predict 
a wide range of possibilities for the large-Q^ region. The upper-left graph of Fig. 4 shows the pre- 
liminary results for the isovector axial form factors using ensembles with 3 different pion masses. 
The data are simultaneously extrapolated in pion mass and (as were the EM form factors). The 
lowest gray band is the result at physical pion mass, and the dashed line is the dipole form with best 
fit to the experimental data {Ma = 1 .026(21 ) GeV [26]). ' The nucleon axial coupling constants are 

^The Ma value used here is obtained from a weighted average of Ma from (quasi)elastic neutrino and antineutrino 
scattering experiments only. Ref. [26] also analyzed weighted values from charged pion electroproduction experiments 
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Figure 5: (upper row) Polarized distribution in a longitudinally polarized proton (left) and neutron (right), 
(lower row) Isovector heUcity distribution of nucleon. 



obtained taking Ga{Q^ = 0) (purple diamonds) which is consistent with other Nf = 2 + 1 results 
(RBC/UKQCD [27] and LHPC [28], statistical eiTorbar only) as functions of m|, as shown in the 
upper-right graph of Fig. 4. The (r^) is obtained from the curvature of the axial form factor, similar 
to the Dirac and Pauli formulation in Eq. 2.2. A dimensionless combination with the lattice nucleon 
mass is plotted as a function of m| in the lower-left graphic of Fig. 4; again, we see consistency 
with other A'y = 2 + 1 results. In the large-m;;: region, we find the lattice products to be consistent 
with the experimental ones, while at lower m^j the deviation becomes more significant. However, 
these points have nijiL < 5; it may be the finite-size effect dominates in this quantity. Further stud- 
ies are needed resolve this issue. Finally, the nucleon axial mass Ma from our calculation and other 
Nf = 2 + 1 ones are displayed in the lower-right of Fig. 4. 

Unlike the nucleon electromagnetic (Dirac and Pauli) form factors, little is known about the 
nucleon axial form factors at large from experiments. With the large-Q^ nucleon axial form 
factors, we can study the polarized distribution in a longitudinally polarized proton (upper left in 
Fig. 5) and neutron (upper right in Fig. 5) for the first time via a similar three-dimensional Fourier 
transformation as Eq. 2.6. The magnitudes of the proton and neutron densities monotonically 
decrease away from the center of the impact plane. Unlike the transverse charge density in the 

and obtained Ma = 1.069(16) GeV. Recent MiniBooNE analysis suggests the Ma could get as high as 1.3 GeV, while 
other recent experiments using neutrino scattering off various targets give a range 1.14-1.26 GeV. Here we took an 
earlier M4 from Ref. [26] as a guideline to compare with our lattice results while the experimental determination of Ma 
remains unclear. 
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q= (GeV=) Q2 (GeV^) 

Figure 6: (left) Experimental data points (seen the legend in the figure) and the monopole (red) and a mixed 
monopole-dipole (dashed) fit to the data. (Right) An example form-factor calculation from the constituent- 
quark model (CQM); different lines (theoretical projections) indicate different quark wavefunctions within 
the same model. Figures taken from Ref [30]. 



unpolarized nucleon, the neutron density does not change sign. Combining with the Dirac and 
Pauli form factors, we can look at the helicity distribution of nucleon in the impact plane [29]; the 
isovector nucleon results are shown in the lower row of Fig. 5. 

2.3 Pion Form Factors 

Recent high-precision pion electroproduction data from Jefferson Lab using L/T-separated 
cross sections provide precise pion form-factor measurements up to 2.45 GeV^ [30]. A summary 
of all the experimental data is shown in the left-hand side of Fig. 6; different-color points indicate 
data from different experiments. See Ref. [30] for more details. The solid red line is the fit using 
the monopole form, and the dashed line is the fit from an 85% monopole -1-15% dipole form to 
the experimental data. Many theoretical models ranging from various quark models to holographic 
QCD provide diverse and disagreeing predictions for the intermediate and high-2^ dependence 
of the pion form factors. Even within the framework of constituent-quark model, differences in 
the treatment of the quark wavefunctions result in very different pion form factors; see the right- 
hand side of Fig. 6 (taken from Ref. [30]). Diverse behaviors are obtained from various theoretical 
calculations. Future precision experiments will challenge these models even more, aiming at 
through 6 GeV^ and higher. 

We apply the same methodology as the nucleon large-Q^ form-factor calculation to the pion 
form factor to understand and make predictions for upcoming experiments. The right-hand side 
of Fig. 7 shows the preliminary results for our dynamical {Nf = 2 + 1) pion form factors with 
580, 875, 1350 MeV pion masses with reaching nearly 7 GeV^ for the highest-mass ensemble. 
The extrapolated form factor at the physical pion mass shows reasonable agreement with JLab 
precision measurements. Future attempts will focus on decreasing the pion masses and exploring 
2^ -dependence of pion form factors for yet higher Q^. 
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Figure 7: (left) Pion form factor with multiple pion masses at 580, 875, 1350 MeV. The experimental points 
are shown as black circles while the lowest gray band is the extrapolation to the physical pion mass using 
our lattice points, (right) bpji from our lattice pion form-factor data. 

We can further use our pion form factor data to study the three-dimensional infinite-momentum- 
frame spatial charge density of the pion [31]: 

P^{b) = r^MbQ)F^{Q^). (2.6) 

JO 271 

Pn and bp]i were first studied using experimental data up to 2.45 GeV^ in Ref. [31], along with 
a few perspectives from perturbative QCD, QCD sum rules, holographic QCD, and the Nambu- 
Jona-Lasinio model. Perturbative QCD predicts the asymptotic behavior for the pion form factor 
at large-2^ region, 

lim F„{Q^) = \6Kas{Q^)fl/Q^. (2.7) 

This leads to a singular central charge density for the pion, which is unusual in nature. 

The right-hand side of Fig. 7 shows bp^ using our calculation of lattice pion form factors. 
Here we use a two-pole form to fit our lattice data. The singularity at = is highly dependent on 
the behavior. However, it is not clear from the current precision which form is better suited 
to describe our data. Further studies are needed to answer this question. 

3. Conclusion and Outlook 

Starting from the earlier exploratory study of large momentum transfer on the lattice for nu- 
cleon Dirac and Pauli form factors, we employed a strategy to overcome noise in the conventional 
method at large Q^. The idea is to include operators that couple to high-momentum and excited 
states and explicitly analyze excited states to get better ground-state signal. We demonstrate that 
one can extend form-factor calculations to higher momentum transfer with reasonable statistical er- 
ror at a fixed lattice spacing. However, systematic error could become significant when one reaches 
a large value of . What we shown for heavier pions on anisotropic lattices can be easily extended 
to isotropic lattice or lighter-pion calculations. 
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Figure 8: Proposed idea for extending the current regime to higher ones 

We demonstrate that our analysis shows consistent (within systematics) results for low-Q^ 
physical quantities, such as {rj 2)' Sa, Ma, comparing with RBC/UKQCD and LHPC's published 
Nf = 2 + 1 results. We calculate larger-Q^ form factors that these calculations lack. With large Q^, 
we can study the various densities of hadrons that depend on a wide range of momentum in form 
factors. In this proceeding, we have shown the results for transverse charge densities in unpolarized 
and polarized nucleons, polarized distributions in longitudinally polarized nucleons, and isovector 
helicity distributions of nucleons. The transverse pion density is also calculated, but we need more 
precision to explore whether there is a singularity at the center of the pion. 

Future work will include smaller source-sink separation to improve the signal; since excited- 
state analysis is included in our methodology, this is not an issue for our ground-state nucleon 
or pion form factors. It would be interesting to study the same form factor using multiple lattice 
spacings to study and reduce systematic errors at large-Q^ due to lattice discretization error. 

To reach even higher Q^, we propose a numerical step-scaling calculation, as shown in Fig. 8), 
On a small volume with very fine lattice spacing, we can easily reach high momentum. By calculat- 
ing the step-scaling function at overlapping momentum points (or interpolating) we can reduce the 
systematic error due to finite-volume or lattice-discretization artifacts. However, generating several 
volumes and tuning for multiple lattice spacing of dynamical Nf = 2 + 1 lattices requires a large 
amount of computational resources; we hope such a proposal will become feasible as petascale 
computing facilities become available in the near future. 
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